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ABSTRACT 

Context. Revised spectroscopic parameters for the HF molecule and a new CN line list in the 2.3 pm region have been recently 
available, allowing a revision of the F content in AGB stars. 

Aims. AGB carbon stars are the only observationally confirmed sources of fluorine. Nowadays there is not a consensus on the relevance 
of AGB stars in its Galactic chemical evolution. The aim of this article is to better constrain the contribution of these stars with a more 
accurate estimate of their fluorine abundances. 

Methods. Using new spectroscopic tools and LTE spectral synthesis, we redetermine fluorine abundances from several HF lines in 
the K-band in a sample of Galactic and extragalactic AGB carbon stars of spectral types N, J and SC spanning a wide range of 
metallicities. 

Results. On average, the new derived fluorine abundances are systematically lower by 0.33 dex with respect to previous determina¬ 
tions. This may derive from a combination of the lower excitation energies of the HF lines and the larger macroturbulence parameters 
used here as well as from the new adopted CN line list. Yet, theoretical nucleosynthesis models in AGB stars agree with the new 
fluorine determinations at solar metallicities. At low metallicities, an agreement between theory and observations can be found by 
handling in a different way the radiative/convective interface at the base of the convective envelope. 

Conclusions. New fluorine spectroscopic measurements agree with theoretical models at low and at solar metallicity. Despite this, 
complementary sources are needed to explain its observed abundance in the solar neighbourhood. 

Key words, stars: AGB and post-AGB; stars: abundances; Physical data and processes: nuclear reactions, nucleosynthesis, abun¬ 
dances 


1. Introduction 


Fluorine has only one stable isotope, 19 F, yet fragile and easily 
destroyed in stellar interiors by proton, neutron and alpha par¬ 
ticle capture reactions. Major fluorine destruction channels are 
19 F(p, cr) 16 0 and l9 F(cr, /?) 22 Ne reactions. Thus, any production 
mechanism has also to enable 19 F to escape from the hot stel¬ 
lar interiors after its production. This makes fluorine a useful 
tracer of the physical conditions prevailing in stellar interiors. 
Fluorine is thought to be produced in several scenarios: (i) Grav¬ 
itational supernovae through the neutrino spallation process, (ii) 
low and intermediate mass (M< 7 M Q ) Asymptotic Giant Branch 
( AGB) star£] during the He-burning Thermal Pulses (TP) and 
the subsequent Third Dredge Up (TDU) episodes, (iii) Wolf- 
Rayet (WR) stars in the hydrostatic He-burning phase and, even, 
(iv) during white dwarfs mergers (see e.g . Woosle^et afllT99CL 
Forestini et al.lfl992t iMevnet & Arnouldll2000fc iLongland et alJ 


201 ill . To date, proof s of their production have been found in 


AGB stars only (e.g. IJorissen et afl 1992; WemeretaL 20051 : 
lOtsuka et al .1120081 lAbia et al]|2010t iLucatello et alJ 2011 ). The 


origin of this element is therefore still rather uncertain and con- 


1 In the following, we define Low Mass Stars (LMS) to be those 
with initial masses between about 0.8 and 2.5 M e which experience 
He ignition under degenerate conditions. Stars more massive than 
this succeed in igniting He gently but they are not sufficiently mas¬ 
sive to ignite C in their core. We named these Intermediate Mass 
Stars (IMS) having initial masses about 2.5-7 M G . 


tributions from the above mentioned sources seem to be nec¬ 
essary when com paring Galactic Chemical Evolution (GC E) 
model predictions (Ren da et al.ll2004l: iKobavashi et al.ll201 lh to 
the evolution i nferred from observations in Galactic field dwar f 
and giant stars (iRecio-Blanco et akll2012t Uonsson et al. 2 014al) . 

To determine the relative role of the above proposed sources 
of fluorine, in particular of AGB stars, accurate determinations 
of F abundances are needed. Fluorine abundances are usually 
determined in cool dwarfs and giants from ro-vibrational HF 
lines at 2.3 pm. This spectral region, unfortunately, is placed 
in a region with a lot of telluric lines that can be difficult 
to be removed. This prevents, in many cases, the use of sev¬ 
eral HF lines to determine the F abundance in a given star. 
The most used one is the R9 A 2.3358 pm HF line, which is 
marginally affected by blends and t e lluric lines (lAbia et al .1 2009 ; 
Ide Lavernv & Recio-Blancoll2013h . IJorissen et al.l (1 1992h firstly 
used this and other secondary HF lines to derive F abundances 
in AGB stars of both oxygen- and carbon-rich type^l showing in¬ 
deed that AGB stars prese nt larg e fluorine enhancements. Most 
recentlv. fAbia et ai . (1200 9, l20l0 ) revisited this pioneering analy¬ 
sis using a new grid of atmosphere models for AGB carbon stars 
and an up-dated molecular and atomic line lists in the 2.3 pm 
region. Yet, they found significant F enhancements but system¬ 
atically lower than those in IJorissen et al.l d 19921) . Nevertheless, 
the new F enhancements are in a global agreement with nucle- 

2 AGB carbon stars are those with surface C/0> 1 (by number). 
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osynthesis calculations of low mass AGB stars with near solar 
metallicity (ICristallo et al.ll20Tlt) . although at lower metallicities 
models predict larger enha ncements than observed (Abi a et alJ 
1201 UlLucatello et all201 ll) . The mentioned abundance analyses 
in AGB stars, however, were performed using HF line excitation 
energies inconsistent with the partition funct i ons adopted in the 
spectral synthesis. In fact, lAbia et all (120091 1201 Ol 201 lh use d 


the Turbospectrum vl0.2 radiative tr ansfer code ( Ple^ 


2012 ), 


1984). 


which adopts partition functions from lSauval & Tatum 
However, the HF line excitation energies used in those works 
were from [Jorissen et alJ (119921) (in turn from Tipping, private 
communication), which are 0.25 eV higher than the correct val¬ 
ues if using the ISauval & Tatuml ( 1984l)’s partition functions. 
As mentioned by Jonss on et alJ (l2014bl) the 0.25 eV different 
adopted excitation energies comes from the fact that the Tipping- 
list uses the dissociation energy of the energy potential, and not, 
like ISauval & Tatuml (119841) . the true energy required for disso¬ 
ciation. The former is higher due to the zero point of the energy 
of the l owest vibrational le vel. The di fference is exactly 0,2 5 eV 
for HF (IZemke et alll99l1) . Recentlv. lJonsson et alJ (l2014bl) pre¬ 
sented a complete and comprehensive calculation of the excita¬ 
tion energies and transition probabilities for the HF molecule 
in the K- and L-bands. The computed HF molecular parame¬ 
ters agree nicely with the new list of experimental molecular pa- 
rameters delivered by the HITRAN database (see iRothman et all 
120131 for details on this database). As a consequence and ac¬ 
cording to the simple curve of growth, F abundances derived 
with consistent partition functions would be lower by an amount 
~ 8Ax, where 9 = 5040/7),// and Ax = 0.25 eV. 

Very recently, a new solar F abundance has been determined 
from the analysis of a sunspot spectrum (iMaiorca et al.ll2014l) . 
The new solar abundance, log e(F)= 4.40 + 0. 2fff is significantly 
lower that the previous commonly adopted value 4.56 + 0.30 al¬ 
though in agreement within the error bars (iAsplund et ali2009l) . 
and it is in a very good agreeme nt with the meteoritic value, 
4.46 ± 0.06 (iLodders et al. 2009). Finally, a new CN line list 
has been provided that might affect the fluorine abundance de¬ 
termination in carbon-rich stars from HF lines in the 2.3 pm 
region (iHedrosa et al.ll20T3l) . As a consequence, in light of the 
new molecular parameters and solar fluorine abundance, a re¬ 
determination of fluorine abundances in AGB stars is needed to 
evaluate their actual ontribution to the Galactic fluorine inven¬ 
tory. This is the aim of the present work. 

In Section 2, we describe the observations and the new spec¬ 
troscopic tools used for the re-determination of the fluorine 
abundances in AGB carbon stars. Section 3 shows our results, 
while we present our conclusions in Section 4. 


2. Analysis and Results 

The stars studie d he re are the same than those analysed in 
lAbia et all (1201(1 and lAbia et al.l (1201 ll) corresponding to Galac¬ 
tic and extragalactic AGB carbon star samples, respectively 
(thereafter Papers I and II). The extragalactic stars belong to 
the Carina dwarf galaxy (2 stars) and the Magellanic Clouds 
(2 stars in the SMC and 1 star in the LMC). Details about the 
observations, data reduction procedures and quality of the final 
spectra can be found in these previously published works. The 
re-determination of fluorine abundances in these stars has been 
made using the spectral synthesis method in LTE with the new 

3 The abundances are given using the usual definition log e(X)= 12+ 
log (X/H), where (X/H) is the abundance of the element X by number 
and log e(H)= 12. 


version (v 14.1) of the Turbospectrum code. Theoretical spectra 
were computed using the same atmospheric parameters (T e ff, 
log g, [Fe/H Q, and microturbulence) as in Papers I & II (see 
those for details), the only difference being the molecular line 
lists used here (however see below about changes in the macro¬ 
turbulence parameter). As noted before, the new HF line exci¬ 
tation energies from [Jonsson et al.l d2014bi) are now consistent 
with the partition functions adopted in the Turbospectrum code. 
The immediate consequence of the systematically 0.25 eV lower 
excitation energies is that, in AGB carbon stars most of the HF 
lines in the 2.3 pm region are saturated (log (W,\lA) > -4.0), 
even for moderate (Solar) fluorine abundances. This makes the 
HF lines more sensitive to the microturbulence parameter (see 
Figure 1) and, as a consequence, the total uncertainty in the de¬ 
rived F abundances would be larger. For instance a change in 
the microturbulence by ±0.2 kins 1 may change the F abun¬ 
dance derived up to +0.15 dex. 

Updates in the molecular line lists in the 2.3 pm region has 
been made. In particular, we have included the contribution from 
the HCN molecule (H I2 CN/H 13 CN) according to the computa¬ 
tions bv lHarris et al.l (120031) . We corrected also the wavelengths 
and intensities of some fFO and OH lines in this spectral re¬ 
gion following the new HITRAN database release. Nevertheless, 
due to the carbon-rich nature of our stars, the contribution to the 
global spectrum of these two later molecules is minimal. This is 
not the case, however, for the HCN molecule, which has some 
contribution by introducing a global extra-absorption (veil) that 
might diminish the spectral continuum up to ~ 2 - 3%. Also, a 
significant change with respect to the line lists used in previous 
works is the updated CN line list in this spectral region. Details 
on how this CN list was computed can be found in lHedrosa et al.l 
( 2013 ). The new CN line list has an impact on the global fitting 
of the observed spectra: theoretical fits are improved for most 
of the stars. For several stars, the C/O ratio corresponding to the 
best fit to the observed spectrum needed to be revised when com¬ 
pared to the derived value in Papers I & II (compare Table 1 in 
these works with Table 2 here). Also, for the stars R Lep and 
SS Vir in Paper I we could not obtain a good fit to their spectra 
and, thus, were discarded from the analysis. On the contrary, we 
find now a good global fit to the spectrum of the extragalactic 
carbon star LMC TRM88 and determine its F abundance (pre¬ 
viously we had just derived an upper limit F abundance for this 
star, see Table 1 in Paperll). The new CN line list also affects the 
HF lines in different ways: while the R9 line still remains almost 
unblended, the R22, R17 and R16 HF lines, which were used in 
lAbia et aD J201Q ) (see their Table 2), are now quite blended with 
CN features (see Figure 1) and have to be considered as sec¬ 
ondary F abundance indicators. On the contrary, the R23, R15 
and R14 HF lines show less blending than before with CN lines 
and complement the R9 line as main fluorine indicators. 

Figure 1 shows examples of synthetic fits compared to ob¬ 
served spectra of the stars TX Psc and UU Aur. It can be 
clearly appreciated the impact of the new excitation energies 
of the HF lines, micro- and macroturbulence parameters and 
that of the new CN line list used in the 2.3 pm region on the 
derivation of the F abundance. 

2.1. The Fluorine Abundance in the Sun and Arcturus 

Solar abundances are used as a zero-point for all the abundance 
studies. Also, Arcturus is a benchmark star in abundance analy- 

4 We use the abundance ratio definition [X/Y]= log (X/Y)*- log 
(X/Y) Q . 
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Fig. 1. Examples of synthetic fits (continuous and dotted lines) to the observed spectra (dots) for typical AGB carbon stars in the region of 
some EIF lines: a) Left. Black line is the best fit to the R9 HF line in TX Psc (log e(F)= 4.7, £ = 2.2 kms 1 and T = 13 kms '); red line is a 
fit assuming a 0.45 dex larger F abundance to show the effect of a F abundance variation similar to the total error estimated on [F/H] (see 
text), while dotted red line shows the corresponding synthetic spectrum adopting the previously wrong lower excitation energy for this 
line, x = 0.48 eV. b) Center. Black line the same as case a; red line is a synthetic spectrum adopting a lower macroturbulence parameter 
T = 11 kms’ 1 , and dotted red line adopting a lower microturbulence, £ = 1.8 kms 1 than in case a, respectively, c) Right. Effect of the new 
CN line list used in this work (black line) compared with the older one (red line) on the fit to the R16 HF line at T ~ 2.2778 fxm in the 
carbon star UU Aur. 


ses of cool giants. To make this analysis as consistent as possible, 
we redetermined the F abundance in these two reference stars 
using the above new molecular parameters and line lists. For the 
Sun, we used the sun spot U mbra Spectra in the region 1.16 to 2.5 
/rm by Wallace e t all (|200 lh . and followed the same procedure as 
in [Majorca et al.l (12014 ) to estimate the effective temperature of 
the atmospheric model that is most compatible with the sunspot 
umbra spectrum. We used weak and moderate intensity OH lines 
at 1.5 //m together with some CO lines at 2.3 gm. Atomic lines 
of the metals present in these spectral regions cannot be used for 
this purpose since most of them are affected by Zeeman split¬ 
ting. This is not the case for the CO and OH molecular lines, 
nor for the HF lines, because the corresponding vibro-rotational 
transitions have no spin nor orbital momentum, which means a 
very weak coupling between the rotational momentum and mag¬ 
netic held. In fact, by using a simple Milne atmosphere model 


and a vertical magnetic field, we checked that these molecular 
lines would be affected by Zeeman splitting only for magnetic 
fields larger than ~ 10 5 6 G, an intensity much larger than that 
observed in the sunspots. We estimated a value T e ff = 3800 
K from the best fit to the CO and OH lines in the sunspot 
spectrum. Thus we used a MARCS model atmosphere with pa¬ 
rameters T e ff/\og g/[Fe/H] = 3800/4.44/0.0 interpolated from a 
grid of atmosphere models JGustafsson et al.l[20081 PI Similarly 


5 We checked this approach by using a more realistic model at¬ 

mosphere fCollados et alJ (19941) for a sunspot which consider the 
effect of magnetic pressure. This sunspot umbra model atmosphere 
is obtained by analyzing Stokes I and V of several spectral lines. 
With the coolest model atmosphere (T~ 3900 K, B= 2750 G at log 
Tgooo = 0.0) deduced by these authors (see their Table 4), which re¬ 
sembles the MARCS model adopted here, we obtain a similar fit to 
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to [Majorca et alJ ( 2014 b we found that the cores of the stronger 
CO and of some OH lines were not so well matched. This dif¬ 
ficulty has already been reported in the analysis of the infrared 
spectra of many K giants with Tpff similar to that we estimated 
in the sunspot (e.g. lRvde et al.ll2002t ITsui3l2008Ll2009h suggest¬ 
ing that the CO fundamental lines cannot be interpreted with a 
photospheric model only. These lines show an excess absorption 
(log W \!A > -4.75) that is probably nonphotospheric in origin, 
indicating that the structure of the hydrostatic atmosphere does 
not describe properly that of the sunsnotiFl Majorca et al.l ( 2014 ’) 
arrived at the same conclusion, nevertheless they estimated a 
T e ff — 4250 K as the best value matching both the OH and 
CO lines in the sunspot spectrum. Nevertheless, we derive log 
e(F)= 4.37 ± 0.04 in the Sun from 10 HF lines, in agreement 
with the value estimated by Maiorcaetal. (l2014l> and with the 
meteoritic abundance (iLodders et al.ll2009)) . The uncertainty on 
this value is mainly determined by the error in the temperature 
estimate of the sunspot spectrum. [Majorca et al.l (1201 4|) quoted 
about 0.25 dex for this error that we also adopt here. In the fol¬ 
lowing we will use our derived F abundance as the reference 
value for the Sun. 

For Arcturus, we used the electronic version of the Infrared 
Atlas Spectru m by HinkleetaL (1995) and the atmosphere pa¬ 
rameters from [Rvde et ali ([2009). We obtain a fluorine abun¬ 
dance of log e(F)= 3.78 ± 0.03 from the analysis of the R7, 
R9 and R12 HF lines. This value is also in e xcell ent agreement 
with t hat obtained by IJonsson et al.l (1201 4al) and IJonsson et al.l 
(l2014bl) . the later from the analysis of some HF lines at 12.2 pm. 


2.2. Fluorine Abundances in AGB Carbon Stars 

Table 2 (second column) shows the fluorine abundances redeter¬ 
mined in the sample of Galactic and extragalactic AGB carbon 
stars of different spectral types. For each star we report the num¬ 
ber of HF lines used (between parenthesis) and the correspond¬ 
ing dispersion when more than a single line was used. Table 2 
also shows (third column) the difference in the F abundance de¬ 
rived with respect to those in Papers I & II. The mean difference 
is -0.33 ± 0.17 dex, i.e. the abundances derived here are on av¬ 
erage systematically lower by this amount. Note that for many 
stars the reduction in the F abundance is different from that ap¬ 
plying just the correction factor ~ 5040/7’,,yyAy (see Table 1). 
The reason is twofold: a) the impact of the new CN line list and, 
b) the use here of larger macroturbulence velocities (T) with re¬ 
spect to those used in Papers I & II. The new CN line list and 
partition functions affects the HF lines in a different manner de¬ 
pending on the actual C/O ratio and the metallicity of the star 
(see Table 1). The former effect depends on the blending of a 
specific HF line with CN lines (see Figure 1) being this, in any 
case, no larger than ~ ±0.20 dex. On the other hand, because the 
use of new molecular line lists (mainly CN), we realised that the 
2.3 pm spectral region is better fitted using a macroturbulence 
parameter in the range 13 - 14 km s' 1 (see Figure 1), instead 
of 11 - 12 kms -1 adopted in Papers I & II. This might increase 

the CO and OH lines and also a very similar average F abundance, 
log e(F)= 4.38 ± 0.04 (see above). This implies that magnetic pres¬ 
sure in a typical sunspot does not affect significantly the derivation 
of the F abundance. 

6 Nevertheless, it is possible to fit reasonably well the cores of these 
molecular lines by extending the atmosphere model up to log Tsooo ~ 
-6. This extension of the photosphere does not affect at all the HF lines, 
which indicates that these lines form into deeper layers where the struc¬ 
ture of our model atmosphere for the sunspot would be more realistic 
(see e.g. lTsuiil2009l for alternative solutions). 



[Fe/H] 


Fig. 2. [F/Fe] vs. [Fe/H] in the stars analysed. Circles data points are: 
black, N-type Galactic carbon stars; red, J-type; green, SC-type and 
blue, extragalactic N-type carbon stars. Dashed lines represent theoret¬ 
ical predictions for a non-rotating 2 M Q TP-AGB star at different TPs 
according to lCristallo et all (2015b). A typical error bar is indicated. 


the fluorine abundance as much as 0.25 dex, compensating par¬ 
tially the systematic decrease of the F abundances because of the 
lowers values. In summary, the F abundances derived here dif¬ 
fer from the previous values by an amount which is different star 
by star depending on the combined effect of the lowers values, 
the higher F parameters and the blending with CN lines, in this 
order of relevance. 

The total uncertainty in the abundance determination can 
be determined from the dependence of the fluorine abundance 
on the stellar atmospheric parameters; this dependence is rather 
similar for all the HF lines, being the uncertainty in T e ff, A F and 
the C/O ratio the main sources of error. A detailed discussion 
on the error sources and their impact in the final F abun¬ 
dance can be found in lAbia et al. ( 2009 ). Also a discussion on 
the uncertainties in the derivation of the stellar metallicity 
[Fe/H ] and v- element abundances can be f ound in Abia et alJ 
d2001il2002h : Ide Lavemv et all (120061 and lAhia et alJ d2008). 
Concerning F, the quadratic addition of all these sources of er¬ 
ror gives a typical uncertainty of ±0.40 dex, together with the 
uncertainty of the continuum position and the dispersion in the F 
abundance when derived from several lines, a conservative total 
error would be ±0.45 dex. This value does not include possible 
systematic errors as those as in the model atmospheres structure 
and/or N-LTE effects. However, the uncertainty in the abundance 
ratio between F and any other element would be lower than this 
value, since some of the above uncertainties cancel out when de¬ 
riving the abundance ratio. Note, however, that the new [F/Fe] 
ratios in many stars (column 5 in Table 2), do not differ signifi¬ 
cantly from the previous ones because the reduction of the solar 
F abundance by 0.19 dex (see previous section). 


3. Discussion 

The main consequence of the redetermined F abundances is that 
they are systematically lower than those in previous determina¬ 
tions, in particular for carbon stars of near solar metallicity. As it 
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Fig. 3. Observed fluorine vs. average s-element enhancements com¬ 
pared with theoretical predictions. Circles symbols are as in Figure 2. 
Solid lines are the predictions for a non-rotating 2 M 0 TP-AGB with 
metallicities (labelled) mimiking those of our stellar sample according 
assuming the standard 13 C-pocket according to lCristallo et al.1 (12015m 
(labelled with F, FRUITY). Dotted lines is the prediction for a non¬ 
rotating 1.3 M q metal-poor model for the same metallicities in the hy¬ 
pothesis of a larger 13 C-pocket (labelled with T, Tail) (see text). Small 
open squares in these lines mark the values predicted at different TPs. 
J-type AGB carbon stars have been excluded because they usually do 
not show s-element enhancements (tAbia & Isemll2000l) . Typical uncer¬ 
tainties are shown. See text for details. 


can be seen in Table 2, normal (N-type) carbon stars show only 
moderate F enhancements ([F/Fe]); the largest enhancements are 
now close to ~ 1 dex but are observed only in one SC-type star 
(WZ Cas) and some of the metal-poor N-type stars in Carina and 
the Magellanic Clouds. On the other hand, J-type carbon stars 
show no F enhancements or, in some cases, even depletions (RY 
Dra and Y Cvn, see Table 2). We confirm, therefore, the results 
of Paper I concerning the [F/Fe] ratios in AGB carbon stars of 
different spectral types. These new enhancements can still be ac¬ 
counted for by current low-mass TP-AGB nucleosynthesis mod¬ 
els (see below). 

We compare the revised F abundances with theoretical mod- 
els computed with the FU NS evolutionary code (Str aniero et alJ 
2006; ICristallo et al.ll2009l) . Those models are available on line 
on the web pages of the FRUITYQ database, which includes 
stars with initial masses 1.3 < M/M 0 < 6.0 and metallic - 
ities -2.15 < [Fe/H] < +0.15 (ICristallo et al] 1201 lL l20l5hh . 
Those models are calculated by coupling the physical evolu¬ 
tion of the structure with a nuclear network including all chem¬ 
ical species, from hydrogen up to Pb-Bi (at the termination 
of the s-process path). Thus, no post process techniques are 
used. In FRUITY models, a thin 13 C pocket forms after each 
TDU episode. The mass extension of the pockets decreases 
along the AGB, thus making the first pockets the most signifi¬ 
cant for the on-going s-process nucleosynthesis. The mass-loss 
rate is calibrated on the observed mass l oss - period relation 
found in Galactic AGB stars (see Stra niero et alJl2006L and ref¬ 
erences therein). In order to properly follow the physical be- 


7 http://fruity.oa-teramo.inaf.it. 


haviour of the most external layers, low temperature C and N 
enhanced molecular opacities are used (ICristallo et al.ll2007l) . In 
fact, if the convective envelope becomes C-rich (i.e. C/0>1), 
C-bearing molecules form largely contributing to the opacity 
(more than O-bearing molecules). Thus, the structure becomes 
more opaque to photons and the external layers expand and cool. 
This implies an increased mass-loss rate. Note that, with re¬ 
spect to previously published FRUITY low mass AGB star yields 
(ICristallo et al.11201 ll) . we find lower F surface abundances due 
to a revision in the opacity treatment of the sub-atmospheric re¬ 


ion (Cristallo et al 


2 015bl), We confirm previous calculations 


(Lugaro et al. 2004; Karakas 2010) showing that the stellar 
mass range for F production peaks around 1.5-2.5 M 0 for all 
Z, showing a small dependence on the metallicity. This im¬ 
plies that very large surface [F/Fe] ratios are predicted for 
decreasing metallicities. Net fluorine yields from IMS stars 
are lower than their low mass counterparts, but indeed have 
a stronger dependence on the metallicity. In fact, in our low Z 
IMS models the 19 F(p, a) ie O and 19 F(a,;?) 22 Ne reactions are 
more efficiently activated, both leading to a net fluorine de¬ 
struction. In any case, however, the IMS contribution to the 
fluorine is negligible, due to the fact t hat these stars ex peri- 
ence a definitely less efficient TDU (see lCristallo et alJ2015bl) . 

Figure 2 shows the [F/Fe] ratios derived in our stars versus 
the stellar metallicity. The different dots refer to the different 
spectral types carbon stars in the Galaxy (see figure caption). 
Blue dots are the extragalactic metal-poor carbon stars belonging 
to the satellite galaxies. From this figure it is evident that J-type 
stars tend to be F-poor objects compared to the normal N-type 
(black dots) carbon stars. Moreover, there is a hint that SC-type 
stars show, on average, larger F enhancements with respect to N- 
type stars. With the present analysis we confirm what we found 
in Paper I in this respect. Note that the origin of SC- and J-type 
stars is still unknown. There is evidence that they are outside the 
classical spectral evolution M-MS-S-C(N) along the AGB phase, 
which is thought to be a consequence of the progressive carbon 
enrichment of the stellar envelope because of TDU episode^ 
Nevertheless, the small number of J- and SC-type stars studied 
prevent us to reach a definite conclusion. We notice from Fig. 
1 that Galactic N-type stars with slightly lower [Fe/H] present 
larger [F/Fe]. This trend is confirmed for the extragalactic AGB 
carbon stars (blue dots), which show definitely larger fluorine 
overabundances, thus confirming the results of Paper II. Dashed 
lines in Figure 2 are the predicted [F/Fe] ratios as a function of 
the metallicity for a typical non-rotating 2 M Q AGB star at differ¬ 
ent TPs (as labelled). It can be seen that observed and predicted 
[F/Fe] ratios agree well within observational errors. Further¬ 
more, the predicted C/O ratios in the envelope at the observed 
metallicity of the stars agree with the corresponding observed 
F enhancement. For instance, in the FRUITY database models 
with [Fe/H]~ 0.0, -1.1, -2.1 (or Z= 0.014,1(T 3 and 10~ 4 , re¬ 
spectively) at the 3 rd TP, the predicted ratios are C/0~ 1.1,2.0 
and 15, respectively. These C/O ratios are compatible with those 
measured in our sample at those metallicities (see Table 2). In 
brief, we confirm here the pre dicted dependence of the F produc¬ 
tion with the metallicity ( e.g. Lug aro et aljl2004t ICristallo et al.l 
l20TIl:lFishlock et akllmfl) . 

A correlation between the F and the ,y-el erne lit overabun¬ 
dance is expected if a large enough 13 C pocket forms after 
each TDU episode. This correlation comes from the 15 N pro- 


8 For a detailed discussion of the propert ies and possible evolu t ionary 
stat us of SC- and J-ty pe carbon stars see Wall erstein & Knannl i ll 9981) 
an d lAbia et all 120031) . 


Article number, page 5 of@ 
















































A&A proofs: manuscript no. fluor 


duction in the radiative 13 C pocket, the site where the s-process 
main component is built up (see Paper I for details). Figure 2 
shows the new [F/Fe] ratios vs. the observed average s-element 
enhancement (Abia et al. 2002; de Laverny et al. 2006; Abia 
et al. 2008). For the solar-like metallicity N-type carbon stars 
(black dots), F and s-element overabundances clearly correlate. 
Such a trend seems to hold for metal-poor extragalactic car¬ 
bon stars (blue dots) as well, with the largest fluorine over¬ 
abundances corresponding to the most s-process rich stars. In 
Figure 3 we compare observations with FRUITY models hav¬ 
ing two different (M,Z) combinations, i.e. initial mass M= 2 
M g with [Fe/H] = 0 and M= 1.3 M 0 with [Fe/H]= -1.67. 
Both curves start from FDU abundances (solid curves). The two 
metalliciteis are representative of galactic and extragalactic car¬ 
bon stars, respectively. At low Z we choose a lower initial mass 
in order to be consistent with the bolometric magnitudes of the 
observed stars (-4.3 <M/,„/ < -5.3). At solar metallicity obser¬ 
vations (excluding most of SC stars) and theory (the two lower 
curves in Figure 3) agree nicely. At low Z (upper continuos 
line in Fig.3), instead, we highlight a clear disagreement, the 
models showing too large 19 F abundances for a fixed s-process 
enrichment (see the discussion in Abia et al. 2011). A simi¬ 
lar apparent deficiency of F abundances for a given s-element 
enhancement has been found also in the Galactic Carbon En¬ 
hanced Metal-Poor stars (CEMP) which are believed to have 
been polluted by mass-transfer from a former AGB star in a 
binary system (e.g. Lucatello et al. 2011). We test whether the 
inclusion of rotation might alleviate this problem by modifying 
the ly F surface distributions in our models (see Piersanti et al. 
2013). The main fluorine production channel starts with the neu¬ 
trons released by the I3 C(cr, «) 16 0 reaction and proceeds through 
the nuclear chain 14 N(n, p) u C(a, y) ls O(p, a) 15 N (a,y ) 19 F (see 
ICristallo et alJ 120141 and references therein). The main effect 
of mixing induced by rotation is to dilute 14 N within the 13 C- 
pocket. Thus, more neutrons are captured by 14 N, feeding the 
nuclear path leading to 19 F production. This has also important 
consequences on the on-going s-process nucleosynthesis, which 
results less efficient depending on the initial rotation velocity 
( 14 N is in fact the major neutron poison of the s-process). Thus, 
rotation leads to slightly larger fluorine surface abundances and 
to a definitely lower s-process efficiency. We conclude, there¬ 
fore, that this physical process cannot be considered as a po¬ 
tential candidate to solve the fluorine discrepancy between the¬ 
ory and observations at low metallicities. A possible way to 
solve such a discrepancy is to obtai n larger s-process enh ance- 
ments for a fixed dredged up mass. ICristallo et al. ( 2015b ) re¬ 
cently demonstrated that a different treatment of the inner bound¬ 
ary of the extra-mixed region at the base of the convective en¬ 
velope during a TDU may lead to appreciable differences in 
the surface s-process distribution. In particular, those authors 
found that allowing the partial mixing to work below the for¬ 
mal Schwarzschild convective boundary down to the layer where 
the convective velocity is 10~ n times lower, a net increase of 
s-process abundances can be obtained. This derives from the 
fact that 13 C pockets are larger than those obtained in standard 
FRUITY models (in which the lower boundary is fixed at 2 pres¬ 
sure scale heights from the formal border of the convective en¬ 
velope). In Figure 3 we report two models with the same com¬ 
bination of mass and metallicity as the FRUITY models already 
discussed, but with the aforedescribed different mixing treatment 
at the base of the convective envelope (dotted lines in Fig. 3). 
At solar-like metallicity we find a slight decrease in the fluo¬ 
rine surface overabundances, this fact however not compromis¬ 
ing the agreeement with observations. At low Z, instead, we no- 



[Fe/H] 


Fig. 4. [F/Fe] vs. [Fe/H] observed in Galactic field dwarfs and gi¬ 
ants and open cluster stars compared with the predicted evolution from 
our GCE model for the solar neighbourhood including only fluorine 
produc tion from AGB stars. B lack line : yield s from ICristallo et all 
d2015ali : Red line : yields from | K arakasl (1 20101) . The observed ratios 
(circle s) are from iRecio- Blanco et all d201 2h (blac k). iLi et all JIoDI) 
(bluet. iMaiorca et afl ~i 2014) (green). I.lonsson et all d2014bl) ('cvanf and 
Nault & Pilachowski (2013) (magenta). A typical error bar is shown. 
See text for details. 


tice a large fluorine reduction for a fixed s-process surface en¬ 
hancement. This leads to a reasonable fit to extra-galactic car¬ 
bon stars. Thus, we conclude that larger L3 C pockets than those 
characterizing FRUITY models are needed to fit fluorine abun¬ 
dances in the metal-poor extra-galactic carbon stars. This fur¬ 
ther strenghten similar conclusions already reached in the s_tudy 
of s-process rich stars belon ging to open clusters (iMaiorca et alJ 
1201 2t iTrin nella et akl l2014l) as wel l as in the is otopic analysis 
of pre-solar SiC grains dLi u et al.112014 l2015li . Note that, in 
any case, SC-type carbon stars marginally agree with theoreti¬ 
cal models. However, the evolutionary status of these objects is 
not well understood. 


3.1. The contribution of AGB stars to the fluorine abundance 
in the Solar neighbourhood 

As mentioned in Section 1, the real contribution of AGB stars 
to the F inventory in the Solar neighbourhood is still a matter 
of debate dJonsson et al.[|2014ah . The present downward revision 
of the F overabundances observed in AGB stars and new model 
predictions re-opens this issue. To shed light on this problem, 
we have computed a simple Galactic Chemical Evolution (GCE) 
model (ICristallo et al.l 12045 a) used recently to study the Solar 
System s-only distribution (i.e. of those isotopes uniquely syn¬ 
thesised by AGB stars). This GCE model accounts for all the 
constraints currently observed in the Solar neighbourhood and 
at the epoch of the Solar System formation (iBoissier & Prantzosl 
119991) . By using this code we have calculated the [F/Fe] vs. 
[Fe/H] evolution including only the F contribution from AGB 
stars. Figure 4 shows the computed evolution (black line) com¬ 
pared to the observed one inferred from Galactic unevolved field 
dwarf and giant stars and open cluster stars (see caption of Fig- 
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ure 4 for the reference^). It is evident that the AGB contribution 
to F is not enough to account for the observed abundance of this 
element in the Solar System: the GCE model predicts a factor 
~ 3 lower [F/Fe] ratio than observed at [Fe/H]~ 0.0. We have 
computed also the evolution using the F yields obtained with the 
STROMLO stellar evolutionary code (see e.g. lKarakasl201(A red 
line). From Figure 4, it is evident that the GCE model produces 
quite different results when using different 19 F yield inputs. The 
reasons for such a discrepancy are manifold. Among the many 
different inputs of the two stellar codes, two main candidates 
can be identified: the treatment of convection and the mas s-loss 
rate. With respect to ours models, those of lKarakasI d2010i ) show 
an increased TDU and hot bottom burning efficiency. While the 
first difference leads to larger l9 F surface enrichments, the sec¬ 
ond implies a more efficient fluorine destruction. However, when 
weighting the contribution with a typical initial mass function, 
the first term becomes dominant. As a consequence, the fluorine 
production in Karakas’s stellar models is larger than our. This 
is particularly evident for stellar models with initail mass M> 3 
M 0 at low metalliciy (Z< 6 x 10 3 ). Another difference arises 
from t he adopted mass-loss rate l aw. The models by Karakas 
(1201 Oh follow the prescriptions of I Vassiliadis & Wood (1 19931) . 
At fixed period, our mass-loss rate is more efficient (see Figure 6 
in Straniero et al. 2006). This leads to a faster erosion of the con¬ 
vective envelope and, thus, to lower yields. Moreover. iKarakasl 
(1201 0b included a metallicity dependency that we do not take 
into account. In any case, using the Karakas (2010) yields for 
F (red line), the computed evolution agrees better with the ob¬ 
served trend at [Fe/H]~ 0.0, but the previous conclusion still 
holds: AGB stars alone cannot account for the current F abun¬ 
dance in the Solar neighbourhood and other sources of F are re¬ 
quired to explain its observed evolution. Note that although us¬ 
ing the Karakas’s yields one might account for the observed 
[F/Fe] ratios at [Fe/H] -0.5 (within the observational un¬ 

certainties, see Fig. 3), at higher metallcities her yields fail to 
get [F/Fe]~ 0.0 as the observations seems to indicate on aver¬ 
age. Additional fluorine abundance determinations in dwarf and 
giant stars covering a wide range of metallicities would help to 
clarify the nature of these additional sources. We notice, further¬ 
more, that current stellar theoretical models reproduce the l<J F 
abundances observed in AGB stars. This further strengthen the 
need of additional fluorine sources. 


4. Summary 

Fluorine abundances have been redetermined in a sample of 
Galactic and extragalactic AGB carbon stars using consistent 
partition functions and spectroscopic parameters of HF lines. A 
redetermination of the F abundance in the Sun and Arcturus has 
been also made, which agree with other recent F abundance de¬ 
terminations in these stars. Concerning AGB carbon stars, the 
new F abundances are systematically lower in average by 0.33 
dex with respect to previous determinations. For a specific car¬ 
bon star the difference depends on a combination between the 
new excitation energies of the HF lines, a higher macroturbu¬ 
lence parameter, and a new CN line list used here. The F abun¬ 
dances found in the near solar metallicity carbon stars agree with 

9 All the [F/Fe] ratios shown in Figure 4 have been re-scaled to the 
new solar F abundance according to Table 2. Note that these ratios 
may not differ significantly from the original ones in the mentioned 
papers, since apart of the lower Solar F abundance used here, one 
should also consider the decrease of the F abundance derived in a 
specific star when using the correct x values. 


the extant nucleosynthesis models in low-mass AGB stars. For 
the metal-poor AGB carbon stars belonging to several satellite 
galaxies, a satisfactory fit can be obtained by using theoreti¬ 
cal models with larger 13 C pockets with respect to our standard 
FRUITY models. Using new F yields from low and intermediate- 
mass AGB stars and a GCE model we evaluate the contribution 
of these stars to the F inventory in the solar neighbourhood and 
conclude that additional sources of F are needed to explain the 
observed evolution of this element. 
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Table 1. Variations (new-old) in the F abundance derived (in dex) from some HF lines due o the combined effect of the news’s and CN molecular 
list for specific atmosphere parameters 


HF lines 


T «// 

logg 

[Fe/H] 

C/O 

R9 

R14 

R15 

2800 

0.0 

0.0 

1.02 

-0.40 

-0.35 

-0.35 

3300 

0.0 

0.0 

1.02 

-0.40 

-0.35 

-0.35 

2800 

0.0 

0.0 

1.10 

-0.40 

-0.35 

-0.40 

2800 

0.0 

-1.0 

1.70 

-0.45 

-0.45 

-0.45 

4000 

3.0 

0.0 

0.58 

-0.33 

-0.30 

-0.25 

5500 

4.0 

0.0 

0.58 

-0.30 

-0.30 

-0.30 


Notes. It is assumed [F/Fe]= 0.0, a macroturbulence T = 13 kms 1 in all the cases, and microturbulence f = 2.2 or 1.5 kms 1 , for carbon- or 
oxygen-rich stars, respectively. 
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Table 2. Fluorine Abundances and Element Ratios 


Star 

log e(F) a 

Alog e(F) to 

C/O 

[Fe/H] 

[F/Fe] 

[<s>/Fe] 

Sun 

4.37+0.04(10) 

0.19 

0.58 

0.00 

0.00 

0.00 

Arcturus 

3.78+0.02(3) 

0.37 

0.20 

-0.53 

-0.07 

-0.01 

N-type 

AQ Sgr 

4.34+0.08 (5) 

0.29 

1.03 

-0.10 

0.08 

_ 

BL Ori 

4.62+0.20 (5) 

0.30 

1.04 

0.00 

0.26 

1.00 

RT Cap 

4.30 

0.36 

1.20 

0.00 

-0.06 

- 

RVCyg 

4.30+0.00 (4) 

0.40 

1.09 

0.00 

-0.06 

0.07 

S Set 

4.37+0.05 (5) 

0.25 

1.05 

-0.10 

0.11 

0.40 

ST Cam 

4.45+0.00 (4) 

0.27 

1.02 

-0.10 

0.19 

0.45 

TU Gem 

4.63+0.24 (5) 

0.29 

1.07 

0.00 

0.27 

- 

TWOph 

4.26+0.06 (3) 

0.35 

1.12 

-0.40 

0.29 

- 

TXPsc 

4.65+0.07 (6) 

0.18 

1.03 

-0.39 

0.69 

1.10 

U Cam 

4.20+0.08 (4) 

0.38 

1.23 

-0.10 

-0.06 

0.55 

UHya 

4.70+0.20 (6) 

0.37 

1.05 

-0.05 

0.39 

1.10 

UU Aur 

4.55+0.14(7) 

0.33 

1.06 

0.00 

0.19 

0.45 

UXDra 

4.56+0.18(5) 

0.29 

1.05 

-0.20 

0.20 

0.65 

V460Cyg 

4.47+0.12(4) 

0.18 

1.07 

-0.05 

0.15 

0.80 

V Aql 

4.26+0.05 (5) 

0.36 

1.05 

-0.05 

-0.05 

0.20 

VY UMa 

4.50+0.03 (5) 

0.20 

1.06 

-0.10 

0.24 

0.90 

WCMa 

4.68+0.13(5) 

0.27 

1.07 

0.20 

0.12 

0.60 

WOri 

4.30+0.00 (3) 

0.12 

1.07 

0.05 

-0.10 

0.25 

X Cnc 

4.32+0.08 (4) 

0.61 

1.03 

-0.30 

0.25 

- 

YHya 

4.40+0.00 (2) 

0.50 

1.18 

-0.10 

0.14 

- 

YTau 

4.36+0.17(4) 

0.21 

1.04 

0.00 

0.00 

0.05 

Z Psc 

4.80+0.06 (5) 

0.10 

1.08 

-0.01 

0.44 

0.97 

SC-type 

CY Cyg 

4.62+0.03 (4) 

0.53 

1.00 

0.10 

0.15 

0.10 

FU Mon 

4.91+0.07 (5) 

0.24 

1.00 

0.00 

0.55 

0.40 

GP Ori 

4.70+0.14(5) 

0.56 

1.01 

0.00 

0.34 

1.25 

RZ Peg 

4.80+0.02 (7) 

0.20 

1.00 

0.00 

0.44 

0.00 

WZCas 

5.31+0.14(7) 

0.49 

1.01 

0.00 

0.95 

0.16 

J-type 

RScl 

4.07+0.05 (4) 

0.11 

1.02 

-0.30 

0.01 

_ 

RYDra 

3.80 

0.68 

1.12 

-0.05 

-0.51 

<0.20 

TLyr 

4.32 

0.35 

1.14 

0.00 

-0.04 

- 

VX And 

4.30+0.00 (2) 

0.05 

1.76 

-0.10 

0.04 

0.2 

YCvn 

< 3.50 

> 0.85 

1.12 

0.00 

< -0.86 

<0.2 

ExtraGalactic N-type 

LMC TRM88 

4.50 

0.35 

1.62 

-0.60 

0.74 

_ 

SMC GM780 

4.10 

- 

3.10 

-0.80 

0.53 

- 

SMC BMB B30 

3.60 

0.00 

1.47 

-1.10 

0.34 

0.90 

Carina ALW6 

3.50 

0.20 

12.0 

-1.70 

0.84 

1.60 

Carina ALW7 

3.80 

0.60 

19.0 

-1.90 

1.34 

1.74 


Note s. [<s>/Fel is the average s-element enhancement according to Abia & Iseml d200flh . lAbia et aTl d2002h . Ide Lavemv et al.1 (2006) and 
lAbiaetalJd2008l) . 

(a> The number in parenthesis indicates the number of HF lines used. (A) Abundance difference in the sense lAbia et aT] d2010h (Galactic stars) 
orjAbiaetaLl ()20nh (ext ragalactic stars) m inus this work. In the cases of the Sun and Arcturus the differences are respect to that derived in 
lAsplund et alJd2009l) and lAbia et alj d2009 ). respectively. 
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